function aaaa_surfemis(lsm, tsu) result(emis)
    use iasing_osse
    implicit none
    ! Arguments
    real(JPR4), intent(in)  :: lsm
    real(JPR4), intent(in)  :: tsu
    integer(JPI4)           :: emis

    if (lsm>0.5) then
        emis = 98
    else 
        if (tsu .lt. T0) then
            emis = 102
        else
            emis = 101
        endif
    endif
end function aaaa_surfemis
